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Abstract. An increase in average air temperature and frequency of rain events is predicted for higher latitudes by 
the end of the 21st century, accompanied by a probable rise in air humidity. We currently lack knowledge on how forest 
trees acclimate to rising air humidity in temperate climates. We analysed the leaf gas exchange, sap flow and growth 
characteristics of hybrid aspen {Populus tremula x P. tremuloides) trees growing at ambient and artificially elevated air 
humidity in an experimental forest plantation situated in the hemiboreal vegetation zone. Humidification manipula- 
tion did not affect the photosynthetic capacity of plants, but did affect stomatal responses: trees growing at elevated 
air humidity had higher stomatal conductance at saturating photosynthetically active radiation (gs sat) and lower 
intrinsic water-use efficiency (IWUE). Reduced stomatal limitation of photosynthesis in trees grown at elevated air 
humidity allowed slightly higher net photosynthesis and relative current-year height increments than in trees at am- 
bient air humidity. Tree responses suggest a mitigating effect of higher air humidity on trees under mild water stress. 
At the same time, trees at higher air humidity demonstrated a reduced sensitivity of IWUE to factors inducing stomatal 
closure and a steeper decline in canopy conductance in response to water deficit, implying higher dehydration risk. 
Despite the mitigating impact of increased air humidity under moderate drought, a future rise in atmospheric humidity 
at high latitudes may be disadvantageous for trees during weather extremes and represents a potential threat in 
hemiboreal forest ecosystems. 
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Introduction 

With rapid increases in global industrial development, fos- 
sil fuel use and changing land-use practices, atmospheric 
CO2 concentration ([CO2]) is expected to double within the 
21st century. This increase will result in global climate 
changes: global mean water vapour concentration, evap- 
oration and precipitation rates, as well as global mean sur- 
face temperature are projected to increase during the 21st 
century (IPCC 2007). These changing climate factors along 
with rising [CO2] affect the physiological performance of 
plants: CO2 assimilation, transpiration, stomatal conduct- 
ance (gs) and ultimately plant growth and productivity. 

The impact of the most common consequences of 
climate change — drought, high temperature and high 
atmospheric vapour pressure deficit (VPD) — on photo- 
synthesis and water use in C3 plants has been quite well 
studied, because the occurrence of extreme tempera- 
tures, soil water deficit and high VPD, as well as their 
interactions, alters the physical properties and yield of 
plants, which are important to agriculture and forestry 
(Fletcher et al. 2007; Guha et al. 2010; Estrada- 
Campuzano et al. 2012; Kuster et al. 2013; Li et al. 2013; 
Sapeta et al. 2013). Considerably less is known of the ef- 
fect of increasing atmospheric humidity on plants. In- 
creases in precipitation are considered very likely at 
high latitudes in the long-term perspective (IPCC 2007). 
Precipitation is predicted to increase in northern Europe, 
especially in winter, and to decrease in southern and cen- 
tral Europe in summer (Rdisdnen et al. 2004). There might 
also be fewer dry days at higher latitudes by the end of 
the 21st century (IPCC 2007). Increasing rainfall fre- 
quency results in higher relative air humidity at local or 
regional scales. 

The leaves of plants grown at high relative humidity 
(RH) have larger stomata, larger stomatal pore aperture 
and length, and significantly lower stomatal density due 
to larger epidermal cells than in plants grown at moder- 
ate RH (Torre et al. 2003; Nejad and Van Meeteren 2005; 
Arve et al. 2013). Therefore, decreasing VPD may lead to 
increased stomatal conductance and to a consequent in- 
crease in transpiration in some plant species grown at 
high RH (Pospisilovd 1996; Fordham et al. 2001; Nejad 
and Van Meeteren 2005). Nevertheless, most findings 
suggest that a decrease in VPD generally leads to de- 
creased steady-state leaf transpiration or sap flux density 
in a wide range of tree species from different habitats 
(Pataki et al. 1998; Meinzer 2003; Bovard et al. 2005; 
Hdlscher et al. 2005). Our previous studies have demon- 
strated decreased sap flux density in response to in- 
creased air humidity in silver birch (6etu/o pendula) and 
hybrid aspen (Populus tremula x P. tremuloides) trees in 
moist summers (Kupper et al. 2011). 



Growing at high RH not only alters stomatal morph- 
ology, but stomatal functioning as well (Fanourakis 
et al. 2010, 2011). It is known that RH is a key environ- 
mental factor mediating changes in stomatal sensitivity 
to CO2 (Talbott et al. 2003). Moreover, RH affects stomatal 
response to water availability and drought. Plants grown 
at high RH are less hydrosensitive than plants grown at 
moderate RH: stomata of high-RH-grown leaves are 
less sensitive to decreases in leaf water potential 
than moderate-RH-grown leaves, and the homogeneity, 
speed and degree of stomatal closure are less in 
high-RH-grown plants (Nejad and Van Meeteren 2005; 
Rezaei Nejad et al. 2006; Rezaei Nejad and Van Meeteren 
2008). Therefore, plants developed under moderate RH 
are able to retain higher water status due to more effi- 
cient stomatal control. Arve et al. (2013) revealed that 
stomata developed under high RH respond to neither 
darkness nor drought, but remain open. Thus, high RH 
may even override the signals given by darkness. The sto- 
mata of plants growing in naturally waterlogged soil are 
also less sensitive to decreasing VPD than those of plants 
growing in well-drained soil (Sellin 2001). 

High RH does not change only the stomatal character- 
istics of plants. Our previous experiments with silver 
birch and hybrid aspen have shown that elevated 
atmospheric RH lowers leaf nutritional status by altering 
nutrient movement via mass flow in soil and lowering 
nutrient transfer through xylem flow into leaves (Tullus 
et al. 20120; Sellin et al. 2013). The changes in leaf 
nutrient content and P:N ratio in turn cause a decline 
in photosynthetic capacity and ultimately changes in 
tree growth rate. 

Experiments on stomatal responses to air humidity 
and plant stress resistance are typically carried out 
in greenhouses or growth chambers with seedlings or 
saplings growing in pots. The objective of the present 
study was to investigate how artificially increased RH 
during leaf development affects the sap flow, stomatal 
responses and photosynthetic parameters of hybrid 
aspen (P. tremula x P. tremuloides) under free-air condi- 
tions. Hybrid aspen is a fast-growing deciduous tree 
species suitable for short-rotation forestry in the rela- 
tively cold climate of northern Europe (Tullus et al. 
2012b). Our aim was to test the following hypotheses, 
(i) Trees grown at higher atmospheric humidity have 
higher stomatal conductance and lower water-use effi- 
ciency (WUE) than control trees, (ii) The photosynthetic 
capacity of leaves developed in humid air is lower be- 
cause of reduced nitrogen uptake due to lower transpira- 
tional flux density, (iii) Plants grown in a more humid 
atmosphere are unable to adjust their WUE quickly be- 
cause of acclimation to lower VPD or possible stomatal 
malfunction. 
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Methods 

study area and sample trees 

Studies were performed on hybrid aspen (P. tremula x 
P. tremuloides) saplings growing in an experimental forest 
plantation at the free-air humidity manipulation (FAHM) 
site, situated at R6l<a village (58°24'N, 27°29'E, 40-48 m 
ASL) in eastern Estonia, representing a hemiboreal vege- 
tation zone. The long-term average annual precipitation 
in the region is 650 mm and the average temperature is 
17.0 °C in July and -6.7 °C in January. In the study year 
(2011) drought conditions prevailed in June and July 
(Fig. 1). The growing season lasts 175-180 days from 
mid-April to October. The soil is a fertile endogenic mollic 
planosol (WRB) with an A-horizon thickness of 27 cm. 
Total nitrogen content is 0.11-0.14 %, C/N ratio is 11.4 
and pH is 5.7-6.3. 

The study site, established on on abandoned agricul- 
tural field in 2006-07, is a fenced area of 2.7 ha contain- 
ing nine circular experimental plots (diameter 14 m) 
planted with hybrid aspen and silver birch (6. pendula) 
and surrounded by a buffer zone. One-year-old micropro- 
pagated hybrid aspen plants were planted in the experi- 
mental area in the autumn of 2006. The stand density in 
the buffer zone is 2500 trees ha"^, and in the experimen- 
tal plots, 10 000 trees ha"^. The computer-operated 
FAHM system, based on an integrated approach of two 
different technologies— a misting technique to atomize/ 
vaporize water and FACE-like technology to mix humidi- 
fied air inside the plots — enables RH of the air to increase 
by up to 18 % over the ambient level during the humidi- 
fication treatment, depending on the wind speed inside 
the experimental stand. The humidificotion is applied 
during daytime 6 days a week throughout the growing 
period if ambient RH is <75 % and mean wind speed is 
<4 m s"^. As a long-term average, RH is 7-8 % greater 
in humidified plots (H treatment) than in control plots 
(C treatment). A detailed description of the FAHM site 




Month of 2011 

Figure 1. Weather data in the growing period of 2011: the dark line 
indicates monthly average air temperature, and the grey bars indi- 
cate monthly precipitation. 



and technical setup is presented in Kupper et al. (2011). 
The treatment began in June 2008; sap flow and gas ex- 
change were measured in the summer months of 2011. 
Soil water potential (^s) was recorded at depths of 15 
and 30 cm with EQ2 equitensiometers (Delta-T Devices, 
Burwell, UK) in eight replications per plot. The daily aver- 
age ^svaried from June to August and was~25 %higher 
in the humidificotion treatment (Fig. 2). The air tempera- 
ture (Ta) and RH were measured 1.5-3.5 m above the 
ground with 2-4 HMP45A sensors (Vaisalo, Helsinki, 
Finland) per plot. Sensor readings were collected every 
1 min and stored as 10-min average values with a data 
logger (DL2e; Delta-T Devices). Air VPD was calculated 
from Ta, saturated vapour pressure and RH. The daily aver- 
age VPD in the humidificotion treatment was 15 % lower 
than the control in the summer of 2011 (Fig. 3). 

Sap flow measurements 

Xylem sap flow in the stems of sample trees was mea- 
sured with FL0W4 sap flow systems (Dynamax Inc., 




Day of year 

Figure 2. Daily average values of soil water potential at a depth of 
15-30 cm in control and humidificotion plots from June to August 
in 2011. Scale bars denote SEM. 
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Figure 3. Hourly average values of air VPD in control and humidifico- 
tion plots from June to August in 2011. Scale bars denote SEM. 



AoB PLANTS www.aobplants.oxfordjournals.org 



© The Authors 2014 



3 



Niglas et al. — 



Hybrid aspen responses to elevated atmospheric humidity 



Houston, TX, USA). Six trees fronn control plots (CI, C2, C4) 
and four trees from humidification plots (HI, H2) were fit- 
ted with sop flow gauges (SGB35-WS) and sampled epi- 
sodically from June to August 2011. Sap flow data were 
recorded every 1 min and stored as 30-min averages. 
Sap flux density (F; g m"^ h"^) was calculated as sap 
flow divided by whole-tree foliage area estimated by 
mean sapwood-to-leaf-area ratios (Huber value, HV) 
measured in nine C (3.08 x 10"^ m^ m"^) and nine H 
(3.09 X lO"'' m^ m"^) trees using destructive sampling. 
Foliage area, measured with a LI 3100C optical area 
meter (LI-COR Biosciences, Lincoln, NE, USA), was on 
average 31 % greater in C trees compared with H trees. 
Whole-tree canopy conductance to water vapour (gc; 
mmol m"^ s"^) was computed from the data of sap flux 
density (mmolm^^s^^) using a simplified Penman- 
Monteith equation (Komatsu et al. 2006; Sellin and 
Lubenets 2010): 

where P is atmospheric pressure (kPa) and VPD is air 
vapour pressure deficit (kPa). 

Gas exchange measurements 

We sampled gasometrically nine trees (mean height 
3.8 m) from C plots and nine trees (mean height 3.3 m) 
from H plots (i.e. three trees per sample plot) for 1 
month, from mid-July to mid-August. Measurements 
were performed on rainless misting-free days on intact 
fully expanded leaves in situ with a portable photosyn- 
thesis system LCpro-i- (ADC BioScientific, Great Amwell, 
UK) at constant air humidity (13 mbar), CO2 concentra- 
tion (Cq = 360 iJimol mor^) and temperature of the leaf 
chamber (25 °C). Leaf-to-air vapour pressure difference 
was relatively similar in the two treatments: on average 
2.12 kPa for C plants and 1.99 kPa for H plants. To gener- 
ate photosynthetic light response curves (A/Q curves), 
four leaves per tree were sampled from the middle part 
of the crown with an instrument equipped with an LED 
light source. The measurements started with photosyn- 
thetically active radiation (PAR) at 1196 (xmol m"^ s"\ 
then decreased stepwise to 9 ixmol m"^ s"^ and increased 
stepwise from 1196 to 1803 jimol m"^ s"^. Intrinsic 
water-use efficiency, expressed as the ratio of net photo- 
synthesis (An) to stomatal conductance to water vapour 
(gs), was determined at two levels of irradiance: at 400- 
600 |xmol m"^ s"^ when IWUE was usually at a maximum 
(IWUEmax) and at light intensities corresponding to full 
sunlight (IWUEsat; 0 > 1400 |xmol m"^ s"^). 

The response of net photosynthesis to varying intercel- 
lular CO2 concentration (C,)— A/Q curves — was also 



determined on intact leaves (four leaves per tree) in situ 
at constant air humidity (13 mbar), temperature of the 
leaf chamber (25 °C) and at saturating irradiance 
(1500 jjLmol m"^ s"^). External CO2 concentration (Cq) was 
supplied in 11 steps, decreasing from 360 to 60 |xmol moP^ 
and then increasing from 450 to 1600 |xmol mor^. In add- 
ition to IWUEmax and IWUEsat calculated from the data of 
A/Q curves, IWUEin (initial IWUE) was determined using 
initial values of the A/Q and A/Q sequences when external 
[CO2] was 360 iJimol moP^. 

Tree growth assessment 

Tree height (H, cm) and stem diameter at 30-cm height 
(D, mm) of all aspen trees growing at three C and 
three H plots were measured before and after the 2011 
growing season. H was measured with a telescopic 
Nedo mEssfix-S measuring rod (Nedo GmbH & Co.KG, 
Dornstetten, Germany) and stem diameter with a LIMIT 
digital caliper (Luna AB, Alingsas, Sweden). Current an- 
nual increment of the trees (AH, AD) was estimated as 
the difference between the two measurements. Relative 
increment (AHrei, ADrei) was expressed as the ratio of AH 
and AD to their respective characteristics at the begin- 
ning of the growing season. The ratio of H : D was defined 
as tree slenderness (S). 

Data analysis 

Statistical data analysis was carried out using Statistica, 
Ver. 7.1 (StatSoft Inc., Tulsa, OK, USA). Repeated- 
measures analysis of variance (ANOVA) was used to com- 
pare the sap flux density (F) and canopy conductance to 
water vapour (gc) between trees from the control and the 
misting treatment. The daily averages of F and gc were 
analysed altogether on 31 days from 1000 to 1700 h 
from 14 June to 7 August 2011 (DOY: 165-176, 197- 
201, 206-219). Linear regression analysis was carried 
out to estimate relationships between F, gc, VPD and ^s- 
The normality of the regression residuals was checked 
using the Shapiro-Wilk test. 

The gasometric data were analysed with Photosyn 
Assistant, Ver. 1.2 software (Dundee Scientific, Dundee, 
UK). The A/Q curves were fitted as a non-rectangular 
hyperbola expressed as a quadratic equation by Prioul 
and Chartier (1977). The initial slope of the curve ex- 
presses the apparent quantum efficiency (</>), whereas 
the X and Yaxes intercepts, respectively, correspond to 
the light compensation point (Ocomp) and apparent dark 
respiration (RJ, and the upper asymptote approximates 
the light-saturated rate of photosynthesis (Amax)- An add- 
itional parameter — convexity (9) — is required to describe 
the rate of bending between the linear increase and the 
maximum value. Sub-stomatal cavity CO2 concentration 
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(Cj) was calculated using the model of von Caemmerer 
and Farquhar (1981). 

The A/Ci curves were analysed according to the bio- 
chemical model proposed by Farquhar et al. (1980), and 
subsequently modified by Harley and Sharkey (1991) 
and Harley et al. (1992). This model enables estimation 
of the CO2 compensation point (F), the maximum rate 
of carboxylation by Rubisco (Vc max), the PAR-saturated 
rate of electron transport Umax) and the rate of triose 
phosphate utilization (Wpu), which indicates the availabil- 
ity of inorganic phosphorus for the Calvin cycle (Sharkey 
1985). The relative stomatal limitation on photosynthesis 



Table 1. Soil water potential (kPa) estimates of the experimental 
plots: 'i's_mean, mean across the growing season; fs 025, lower 
quartile; ¥5 075, upper quartile. 



Plot 


^S_nnean 






CI 


-197 


-240 


-152 


C2 


-191 


-217 


-177 


C4 


-196 


-217 


-185 


HI 


-56 


-76 


-33 


H2 


-194 


-221 


-175 


H4 


-124 


-151 


-97 



(Ls), an estimate of the proportion of the reduction in 
photosynthesis attributable to CO2 diffusion between 
the atmosphere and intercellular space, was calculated 
from the A/Q curves as follows (Farquhar and Sharkey 
1982; Tissue et al. 2005; Huang et al. 2008): 




where An is the net photosynthetic rate at normal Cq 
(360 p,mol mor^) and Aq is the photosynthetic rate 
when Ci ( = 360 |xmol mor^) equals Cq. Under these con- 
ditions, Ao is the rate of photosynthesis that would occur 
if there were no diffusive limitation to CO2 transfer 
through stomatal pores. The effect of humidification on 
gas exchange parameters was analysed by applying a 
nested analysis of variance with fixed factors of 'Treat- 
ment', 'Experimental plot' and 'Soil water potential' (a 
continuous variable), the second nested in the first. As 
plant physiological traits were more strongly related to 
the soil water potential measured at 30-cm depth (^30), 
we used this parameter as an index of soil water status. 
Because of drought development during the measure- 
ment period, we divided the datasets of both treatments 
into two groups according to (< -204 kPa for drier 




10 11 12 13 14 15 16 17 10 11 12 13 14 15 16 17 



Time of day Time of day 

Figure 4. Daily average values of sap flux density (F) and canopy conductance to water vapour (gd in control and humidification plots during 
mist fumigation from June 14 to June 25 (DOY: 165-176; A and C) and July 15 to August 7 (DOY: 197-219; B and D), 2011. Scale bars denote SEM. 
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soil and > - 204 kPa for nnoister soil in C plots; < - 163 
and >163 kPa in H plots, respectively) and analysed gas 
exchange data also separately for these conditions. 

Student's t-test was applied to estimate the treatment 
effect on the growth characteristics of individual trees 
across all experimental plots. Analysis of variance models 



o 

CJ 

a: 

o 
c 
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O 
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SO 
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Figure 5. Variation in daily average canopy conductance to water 
vapour (gc) depending on atmospheric VPD (A) and bulk soil water 
potential ata depth of 30 cm ('i'Bo; B).The numbers by the regression 
lines indicate the respective slopes. 



were used to study the effects of 'Treatment' and 'Experi- 
mental plot' (nested in treatment) or 'Treatment' and 'Soil 
water potential' as a continuous covariate on the growth 
characteristics. Means and upper and lower quartiles of 
daily average soil water potentials (^s mean, ^s_q25, 
^s_Q75) across the growing season were used as covari- 
ates in separate models (Table 1). When exploring the 
variance of total and relative growth increment in 2011, 
tree size (H or D at the end of the previous growing sea- 
son) was included as a covariate. Type IV sums of squares 
were used in the calculations; post hoc mean compari- 
sons were conducted using Tukey's HSD test. 

Results 

Sap flux density and canopy conductance 

Although canopy conductance (gc) was significantly 
higher (22 %; P < 0.05) under humidification across the 
whole study period, the difference between the treat- 
ments was statistically insignificant for days 165-176 
(Fig. 4C) when the soil water potential did not differ be- 
tween the C and H plots (Fig. 2). Also the sop flux density 
in the H treatment was on overage 13 % higher than in 
the C treatment, although the difference was not signifi- 
cant (Fig. 4A and B). gc decreased with increasing VPD 
(P< 0.001) in both the C and H plots; the response pat- 
terns were completely coincident and the slopes of the 
respective regression lines did not differ between the 
treatments (Fig. 5A). gc also decreased with decreasing 
"^s (P < 0.001), while the treatments demonstrated con- 
trasting sensitivities (dgc/d^s) to developing soil water 
deficit — the corresponding slopes were 0.94 and 3.01 
for control and humidified trees, respectively (Fig. 5B). 



Table 2. Leaf gas exchange characteristics of hybrid aspen growing under control and humidification treatment. Each value is the mean + SE; 
the means are compared with Tukey's test. NS, not statistically significant. 



Parameter 


Treatment 
Control + SE 


Humidification + SE 


Significance level (P) 


An (jjLmol m"^ s"^) 


10.5 


+ 0.4 


11.3 + 


0.4 


NS 


Qs (mol m"^ s"^) 


0.22 


+ 0.01 


0.28 + 


0.01 


0.002 


Qs sat (mol m"^ s"^) 


0.19 


+ 0.02 


0.25 + 


0.02 


0.040 


IWUEin (ixmol mol"^) 


53.4 


+ 1.5 


45.1 + 


1.5 


<0.001 


IWUEmax (jj-molmol"^) 


62.3 


+ 2.29 


54.4 + 


2.19 


0.015 


IWUEsot (M-mol mol"^) 


56.8 


+ 1.88 


48.81 + 


2.19 


0.008 


Ls (%) 


41.3 


+ 1.01 


37.7 + 


1.01 


0.016 


Amax (fJLmol m"^ s"^) 


12.9 


+ 0.6 


12.9 + 


0.6 


NS 


Vc max (iJLmol m"^ s"^) 


56.5 


+ 2.9 


59.9 + 


2.7 


NS 


Jmax (fj-mol m"^ s"^) 


173 


+ 10 


196 + 


11 


NS 
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Leaf gas exchange 

Average net photosynthesis (An) tended to be slightly 
greater in trees growing at elevated atmospheric humid- 
ity than those grown at ambient RH, although the treat- 
ment means did not differ statistically throughout the 
experiment (Table 2). Analysis of variance revealed that 
the humidity treatment affected stomatal response, but 
not leaf photosynthetic traits (Tables 2 and 3). Specifically, 
there were significant differences in means of stomatal 
conductance to water vapour measured at saturating 



PAR (Qs sat) and IWUE between the treatments: gs sat was 
32 % higher and IWUEin 16 % lower in the H treatment 
than in C trees (Table 2). 

The data analysis revealed that soil water availability 
affected the gas exchange parameters differently 
within the treatments. and gs sat in the H treatment 
were significantly greater under moist soil conditions 
(12.45 |xmol m"^ s"^ and 0.300 mol m"^ s"\ respectively) 
than under drier conditions (9.78 |xmol m"^ s"^ and 
0.192 mol m"^ s"^, respectively; Fig. 6A and B). Initial 



Table 3. Effects of treatment, plot and soil water potential at a depth of 30 cm {Wsol on gas exchange characteristics. NS, not statistically 
significant. 



Characteristic 


Factor 


Significance level (P) 


An (jjLmol m~^ s"^) 


Treatment 


NS 




Plot (nested in treatment) 


< 0.001 






0.013 


gs (mol m"^ s"^) 


Treatment 


NS 




Plot (nested in treatment) 


0.002 




^30 


NS 


gs sat (molm"^s"^) 


Treatment 


0.027 




Plot (nested in treatment) 


0.022 




*'30 


NS 


IWUEin (M-mol mol"^) 


Treatment 


0.013 




Plot (nested in treatment) 


0.002 




*'30 


NS 


IWUEmax (M-mol mol"^) 


Treatment 


< 0.001 




Plot (nested in treatment) 


0.005 




^30 


0.003 


IWUEsat (M-mol mol"^) 


Treatment 


< 0.001 




Plot (nested in treatment) 


0.012 




^30 


0.007 


^max 


Treatment 


NS 




Plot (nested in treatment) 


NS 




*'30 


NS 


Vc max (|Ji-mol m"^ s"^) 


Treatment 


NS 




Plot (nested in treatment) 


NS 




*'30 


0.045 


Jmax (fJLmol ITn"^ S"^) 


Treatment 


NS 




Plot (nested in treatment) 


NS 




^30 


0.007 


Ls (%) 


Treatment 


NS 




Plot (nested in treatment) 


< 0.001 




^30 


0.015 
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Figure 6. Means of net photosynthesis (Ar,; A), stomotal conduct- 
ance to water vapour at saturating PAR (g^ sat) B) and intrinsic water- 
use efficiency (IWUEir,; C) of control (closed circles) and humidified 
trees (open squares) depending on soil water status. Values are 
means + SE; different letters denote statistically significant 
(P < 0.05) differences. 



intrinsic water-use efficiency also differed with respect to 
soil conditions: it was lower under moist soil conditions 
than under drier conditions (44.2 and 49.3 p,mol mor^, 
respectively). There were no differences in An, Qs sat and 
IWUEin with respect to soil moisture in C plots (Fig. 6). It 
is important to notice that gs, sat. IWUEin and Amax did 
not depend on ^so (as a continuous variable) across the 
whole dataset (Table 3). 

Photosynthesis was strongly associated with gs in both 
treatments: = 0.84, P < 0.001 in C plots and = 0.79, 
P < 0.001 in H plots. There was an inverse linear relation- 
ship between IWUEm and Q/Ca (R^ = 0.73, P< 0.001), 



while the slopes of the corresponding regressions did 
not differ between the treatments (P > 0.05). An in- 
creased with rising CO2 concentration (Cq), with signifi- 
cantly (P < 0.001) steeper response in the H treatment 

= 53.5; R^ = 0.82, P < 0.001) than in the control {p = 
45.5; R^ = 0.74, P < 0.001). There were no differences in 
the An = f(Ci) slopes between the treatments, gs de- 
creased with increasing Cq, but the responses did not dif- 
fer between the treatments. As a consequence, IWUEjn 
rose with Ca, the slope for control trees was greater 
than that for humidified trees (P < 0.001; Fig. 7). 

The maximum rate of carboxylation by Rubisco (Vc max) 
and the maximum rate of electron transport Umax) did 
not differ between the treatments across the whole data- 
set; ANOVA revealed only an effect of soil water status on 
these parameters (Table 3). When the data were analysed 
separately in two groups (moist versus dry soil condi- 
tions), significant differences between the means of 
max and Jmax become evident only for the humidifica- 
tion treatment — both parameters were higher (P < 0.001 
for both parameters) in moist soil. No variation with 
the soil conditions was detected in Vc max and Jmax in 
the control trees (Fig. 8A and B). Regression analysis re- 
vealed a positive relationship between Vc max (R^ = 0.168, 
P < 0.001) and Jmax (R^ = 0.151, P < 0.01) and ^^30, as 
well as between Vc max and Jmax (R^ = 0.85, P < 0.001) 
across both treatments. 

The mean values of relative stomatal limitation 
of photosynthesis (Ls) were lower in trees grown at ele- 
vated RH than in C trees— 37.9 and 41.3 %, respectively 
(P< 0.05; Table 2), although ANOVA did not establish 
any significant effect of the treatment (Table 3). Net 
photosynthesis was negatively correlated with Ls in con- 
trol trees (R^ = 0.15, P = 0.03), but the relationship lacked 
in the humidification treatment (P = 0.23). We found no 
differences in Ls with respect to soil water status in any 
treatment separately (Fig. 8C). 



Impact on growth rote 

Saplings of hybrid aspen growing in H plots were signifi- 
cantly shorter and had narrower stems (Table 4), regard- 
less of whether sample plot or soil water potential was 
included as confounding factors in the models (Table 5). 
The absolute and relative growth increments in 2011 
were either unaffected by treatment or significantly 
greater in H plots (Tables 4 and 5). This was more pro- 
nounced when ^s_mean Of ^s_Q25 wos used OS 0 covoriote 
in ANOVA models, although using ^s_q25 yielded slightly 
better approximations than the two other soil water po- 
tential estimates (Table 5). Slenderness (S) was unaffected 
by treatment, but varied significantly among the experi- 
mental plots. 
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Discussion 

Effects on sap flow and gas exchange 

The sap flux density in hybrid aspen trees changed con- 
siderably compared with that at the same experimental 
site in previous summers (Kupper et al. 2011; Tullus 
et al. 2012o). F in control plots was significantly greater 
than that in the humidificotion treatment during the 
rainy summer of 2009 (Kupper eto/. 2011). The some ten- 
dency (not significant) was observed in the drier summer 
of 2010 (Tullus et al. 2012a). However, our current results 
demonstrate higher (although not statistically signifi- 
cant) sap flux densities in trees growing in the humidifico- 
tion plots (Fig. 4). The canopy conductance to water 
vapour in the H treatment was greater (P < 0.05) across 
the whole study period. This discrepancy is attributable 
to relatively low soil water potential in the control treat- 
ment owing to the very dry summer: did not rise over 
-80 kPa in July. The total precipitation during the 2011 
growing season (May-October) was 261 mm (Fig. 1), 
42 % less than the overage of the three previous years 
(452 mm). As the air humidity manipulation did not affect 
sopwood-to-leaf-oreo ratio (HV), the ~30 % greater leaf 
area of control trees was responsible for greater tron- 
spirotionol water losses causing faster depletion of soil 
water reserves and a greater decline in despite lower 
overall sop flux densities in C plots (Fig. 4A and B). 

The response of canopy conductance to changes in VPD 
did not vary between the treatments, suggesting that 
stomotol sensitivity to atmospheric evaporative demand 
was not affected by the experimental manipulation. In 
contrast to that, gc decreased much faster in response 



to falling in the H treatment than in the control 
(Fig. 5). The differential response of gc to decreasing soil 
water availability is probably mediated by plant hydraulic 
conductance (Cohen and Naor 2002; Domec et al. 2009). 
Hydraulic measurements performed on aspen trees in 
2010 revealed that both soil-to-branch and leaf hydraulic 
conductances expressed per unit leaf area were smaller 
in humidified trees, although growing in moister soil 
(A. Sellin, unpubl. res.). Under conditions of soil water 
deficit the lower hydraulic capacity probably becomes 
a crucial factor for H trees, limiting leaf water supply 
and inducing a steep decline in canopy conductance. Fur- 
thermore, on experiment with silver birch revealed that a 
rapid water deficit in H plants led to a faster decrease in 
hydraulic conductance — responsible for liquid water sup- 
ply—compared with the decrease in gs, which limits 
water losses, and exposed plants to a greater risk of de- 
hydration (Sellin et al. 2014). 

Lowering and gs and increasing IWUE ore typical re- 
sponses to water stress in plant species with a drought 
avoidance strategy. When plants encounter a soil water 
deficit, obscisic acid (ABA) is synthesized in the roots 
and translocated to the leaf through the transpiration 
stream (Assmann and Shimozaki 1999); higher concen- 
trations of ABA in leaves drive mechanisms leading to a 
decrease in gs and on increase in WUE (Liu et al. 2005). 
Pontin et al. (2013) propose that ABA promotes stomotol 
closure in two ways— via its widely known biochemical ef- 
fect on guard cells and via on indirect hydraulic effect 
through a decrease in leaf hydraulic conductance. Main- 
taining stable gas exchange attributes during drought de- 
velopment means that o plant either possesses a drought 
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Figure 8. Means of maximum rate of carboxylation by Rubisco 
(Vc max! A), PAR-scturoted rate of electron transport (JmaxJ B) and 
relative stomatal limitation to photosynthesis (Ls; C) of control 
(closed circles) and humidified trees (open squares) depending on 
soil water status. Values are means + SE; different letters denote 
statistically significant (P< 0.05) differences. 

tolerance strategy or lacks adaptations with respect to 
drought. Populus tremula, one of the parent species, is 
known to act as a drought avoider (Possen et al. 2011). 
In our case, the lack of variation in gas exchange charac- 
teristics in C plots with respect to soil water status (Fig. 6) 
and modification of gas exchange in H plots by altering 
water availability can be explained by the lower values 
in the C treatment (i.e. long-term effects mediated 
by ABA). 



Regardless of the difference in between the treat- 
ments, there was still an effect of the humidity manipula- 
tion on Qs sat and IWUE (Table 3). In fact, the differences in 
leaf gas exchange between the C and H plots are attrib- 
utable to the combined effects of soil water availability 
and increased atmospheric humidity. As such, our first 
hypothesis is supported by the experiment: growing at 
higher RH increases stomatal conductance in trees 
while lowering photosynthetic WUE (Table 2), while the 
effect is largely mediated by changes in soil water status. 

Soil water potential influenced both Vc max and J^ax in 
hybrid aspen (Table 3). Grassi et al (2005) found a positive 
relationship between Vc max and in oak and ash trees 
during summer, as in hybrid aspen in this study. A simu- 
lation by Grassi and Magnani (2005) indicated that 
30-40% of the biochemical limitation could be attribu- 
ted to a reduction in leaf nitrogen content during 
droughty summers. Plant photosynthetic capacity and 
leaf N content expressed per leaf area are positively 
correlated (Grassi et al. 2005; van de Weg et al. 2012). 
Previous work at the FAHM experimental site has shown 
that rising RH reduces the water flux in silver birch 
(Kupper et al. 2011) and alters the nutritional status of 
leaves, leading to a decline in photosynthetic capacity 
(Sellin et al. 2013). In this experiment, humidification 
increased rather than decreased water flux through 
the trees (Fig. 4), which explains why the biochemical 
capacity of photosynthesis was unaffected by the 
manipulation, but by changes in (Table 3). The second 
hypothesis concerning reduced photosynthetic capacity 
of leaves due to the expectedly smaller N uptake in 
humidified trees remained unconfirmed. The absence 
of on impact of air humidity manipulation on photosyn- 
thetic machinery of hybrid aspen was also supported by 
chlorophyll fluorescence measurements performed in 
the droughty summer of 2012 (A. Niglas, unpubl. res.). 

High RH can affect stomatal sensitivity by changing 
stomatal morphology: plants grown at higher RH hove 
larger stomata that close to a lesser extent when leaves 
dry (Giday et al. 2013). In addition, long-term acclimation 
to high RH during growth increases heterogeneity in sto- 
matal response characteristics to short-term exposure 
to stomatal closure-inducing factors (Nejad and Van 
Meeteren 2005). Experiments showing differences in sto- 
matal sensitivity and morphology between plants grown 
at high and low RH have been carried out under stable/ 
controlled environmental conditions. The conditions be- 
fore a leaf is fully expanded are important determinants 
on whether stomatal closure capacity is affected by leaf 
dehydration and RH. Moreover, the degree of stomatal 
adaptation in expanding leaves depends on the duration 
and timing of the exposure to high RH (Fanourakis et al. 
2011). The present study was performed under field 
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Table 4. Comparison (t-test) of mean ( + SE) growth parameters of individual trees growing in humidified (H) and control (C) plots. NS, not 
statistically significant. 



Growth choroctcristic 


r 
L 




u 
n 


r-siai 


n 
r 


Height (cm) 


454 


+ 5.8 


428 + 5.8 


3.06 


0.002 


Diameter of stem at 30-cm height (mm) 


32.6 


+ 0.7 


29.5 + 0.6 


3.39 


<0.001 


Height increment of the current year (cm year" ■') 


111 


+ 2.3 


117 + 2.7 


-1.70 


NS 


Relative height increment of the current year 


0.34 


+ 0.01 


0.39 + 0.01 


-3.46 


<0.001 


Diameter increment of the current year (mm year"^) 


8.9 


+ 0.3 


8.5 + 0.2 


1.14 


NS 


Relative diameter increment of the current year 


0.38 


+ 0.01 


0.41 + 0.01 


-1.65 


NS 


Slenderness (height-to-diameter ratio) 


14.4 


+ 0.2 


14.8 + 0.2 


-1.97 


0.051 



Table 5. Results from ANOVA models describing the effect of tree size (i.e. the value of the respective parameter before the start of the growing 
season, Q-i), treatment (T), plot and soil water potential (iJ's mean. 'i's_Q25. ^s_q75) on the growth parameters. 



Factors 




Response variables 












H 


D 


AH 




AD 


ADrel 


s 


Model 1 


















Q-i 


P 






0.824 


<0.001 


< 0.001 


0.766 




T 


P 


< 0.001 


<0.001 


0.074 


0.247 


0.069 


0.086 


0.051 


Plot (7) 


P 


< 0.001 


<0.001 


<0.001 


<0.001 


< 0.001 


< 0.001 


0.015 




Adj. 


0.29 


0.13 


0.09 


0.37 


0.37 


0.16 


0.03 


Model 2 


















Q-i 


P 






0.802 


<0.001 


< 0.001 


0.006 




T 


P 


< 0.001 


<0.001 


0.019 


0.013 


0.041 


0.017 


0.668 


^Smean 


P 


< 0.001 


0.157 


<0.001 


<0.001 


< 0.001 


< 0.001 


0.005 




Adj. 


0.08 


0.03 


0.09 


0.36 


0.32 


0.09 


0.03 


Model 3 


















C,-i 


P 






0.880 


<0.001 


< 0.001 


0.004 




T 


P 


< 0.001 


<0.001 


0.013 


<0.001 


0.049 


0.018 


0.800 




P 


< 0.001 


0.081 


<0.001 


<0.001 


< 0.001 


< 0.001 


0.010 




Adj. R^ 


0.09 


0.04 


0.10 


0.37 


0.32 


0.09 


0.03 


Model 4 


















Ct-i 


P 






0.684 


<0.001 


< 0.001 


0.008 




T 


P 


< 0.001 


0.002 


0.083 


0.077 


0.075 


0.052 


0.622 


fs_Q75 


P 


< 0.001 


0.351 


<0.001 


<0.001 


< 0.001 


< 0.001 


0.002 




Adj. R^ 


0.07 


0.03 


0.07 


0.35 


0. 32 


0.08 


0.04 



conditions with natural diurnal fluctuations of RH; misting 
was applied when ambient RH was <75 % and could be 
increased to as much as 18 % (versus 60 and 95 % of RH 
in Fanourakis et al. 2011). Our data suggest that stomatat 
sensitivity to atmospheric VPD remained unaffected 
in saplings of hybrid aspen. Although we did not explore 



stomatal dimensions, we presume that differences 
in morphology and the putative morphological effect 
on stomatal sensitivity were rather minor as our trees 
grew in natura, in both diumally and seasonally variable 
environments, under conditions requiring flexible stomatal 
adjustment. 
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The findings that stomatal conductance decreases and 
photosynthesis increases with rising external CO2 level 
are well-known phenomena (reviewed by Araujo et al. 
2011). A steeper An response to Cq in H plots is attribut- 
able to higher stomatal conductance (evidenced by gs 
max) due to leaf development under lower VPD (Table 2). 
This Is indirectly confirmed also by the negative correl- 
ation between An and Ls In the C treatment. 

Initial intrinsic water-use efficiency responded more 
sensitively to Cq in C plots than in H plots (Fig. 7), testi- 
fying once more to the effect of elevated atmospheric 
humidity on leaf gas exchange. Because high IWUE Is 
advantageous to plants under drought conditions, a 
slower response in hIgh-RH-grown trees to changing 
ambient conditions may be disadvantageous In the 
case of abrupt climatic fluctuations becoming more fre- 
quent in the future (Easteding et al. 2000); these plants 
are not able to adjust their water use as quickly as plants 
grown in drier air and experience greater water loss. 
Thus, our results support the third hypothesis on the 
capacity of plants to modify WUE under changing envir- 
onmental conditions, albeit not directly tested with 
respect to air humidity. 

Consequences on tree growth 

The above-ground growth response of aspen trees to 
humldificotion in 2011 demonstrated some trends in con- 
trast to those observed in previous years (Tullus et al. 
2012a). However, the positive effect of humldificotion 
was detectable only In current-year height increments, 
while overall dimensions remained smaller In H plots, 
where hybrid aspen trees had grown slower than in 
C plots in the two previous experimental years (Tullus 
et al. 20120). The Inverse growth response Is also attribut- 
able to dry weather conditions prevailing In summer 
2011. Generally or experimental plot was the 
more significant factor influencing tree growth response 
than humidity manipulation. One must take into account 
that the two factors— ^fs and treatment— are partly in- 
terrelated, as transplratlonal flux through trees was 
lower in H plots (Fig. 4A; see also Kupper et al. 2011 and 
Tullus et al. 2012o) and more water was retained In the 
soil (Fig. 2). However, the humidity manipulation also 
had an impact on growth when considering the effect 
of thus, the humldificotion effect on tree growth 
was clearly not due solely to altered soil water availability. 
In average or rainy years, when soil water does not limit 
growth, lowered transpiration hinders nutrient uptake by 
trees in H plots (Tullus et al. 2012o), especially for nutri- 
ents migrating to the roots with mass flow in soil. Under 
these conditions increased atmospheric humidity does 
not improve the growth rate of hybrid aspen. Sellln et al. 
(2013) also showed that humldificotion treatment lowers 



the photosynthetic capacity and growth rate of silver 
birch In moist summers. In dry years, when soil water 
availability limits growth, the impact of this mechanism 
Is obviously Irrelevant. 

Conclusions 

The current study demonstrates that higher air humidity 
mitigates the effect of low soil water availability on 
broadleaved trees during dry years by reducing stomatal 
limitation to photosynthesis, allowing higher net photo- 
synthetic rotes and supporting higher growth rates 
(relative height growth). At the same time, higher RH In- 
creases the sensitivity of canopy conductance to water 
deficit and reduces the responsiveness of IWUE to factors 
Inducing stomatal closure. The present and our earlier re- 
sults (Tullus et al. 2012o; Sellln et al. 2013, 2014) imply 
that 0 future rise In atmospheric humidity at high lati- 
tudes may be disadvantageous In evenly rolny/humld 
years and expose trees to a higher dehydration risk during 
weather extremes, although mitigating the Impact of soil 
water deficit under moderate drought. 
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